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An efficient broadband sum frequency generation (SFG) technique using the two collinear optical parametric 
processes w% = ui\ + u>2 and 0J4 = oji + u)s is proposed. The technique uses chirped quasi-phase-matched 
gratings, which, in the undepleted pump approximation, make SFG analogous to adiabatic population transfer 
in three-state systems with crossing energies in quantum physics. If the local modulation period first makes 
the phase match occur for UI3 and then for U4 SFG processes then the energy is converted adiabatically to the 
UJ4 field. Efficient SFG of the LJ4 field is also possible by the opposite direction of the local modulation sweep; 
then transient SFG of the field is strongly reduced. Most of these features remain valid in the nonlinear 
regime of depleted pump. © 2011 Optical Society of America 
OCIS codes: 190.4223,130.4310,190.7220,270.1670. 



Introduction. Recent advances in quasi-phase- 
matching (QPM) techniques [1, 2] have drawn analo- 
gies between optical parametric processes and two- and 
three-state quantum systems [3-5]. By using an anal- 
ogy to stimulated Raman adiabatic passage (STIRAP) 
in atomic physics [6-11] Longhi proposed [3] a scheme 
in which the fundamental frequency field is directly con- 
verted into the third harmonic without a transient gener- 
ation of the second harmonic. This proposal requires the 
simultaneous phase matching of second harmonic gen- 
eration (lu + lu = 2u>) and sum frequency generation 
(lu + 2lu = 3uj); because this condition can be fulfilled 
only for a specific frequency this technique is not broad- 
band. Suchowski et al. [4, 5] used an aperiodically poled 
QPM crystal to achieve both high efficiency and large 
bandwidth in sum frequency generation (SFG) in the 
undepleted pump approximation using ideas from rapid 
adiabatic passage in quantum physics [8, 12]. 

In this Letter, we make use of the analogy between 
coherent population transfer in three-state quantum sys- 
tems and the two simultaneous collinear second-order 
parametric processes cu 3 = u)\ + ui2 and LU4 = lui + uj 3 ~ 
2lu\ + LJ2 to design a potentially highly efficient broad- 
band SFG technique. To this end, we use linearly chirped 
QPM gratings [13-15] , which provide the analogy to level 
crossings in atomic systems [16-18]. 

The two simultaneous SFG processes lo 3 = oj\ + 1^2 
and lu 4 = oj\ + lu 3 , for a QPM crystal with susceptibility 
and local modulation period A(z) are described by 
the set of nonlinear differential equations [1,2] 
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where z is the position along the propagation axis, c 
is the speed of light in vacuum, and Ej, LOj and rij 
are the electric field, the frequency and the refractive 



index of the j-th laser beam, respectively. Here f2j = 
X^Wj/Acrij (j = 1,2,3,4) are the coupling coefficients, 
while Ai ~rii/ (cuji) + n 2 / {0x12) ~ n 3/ (CW3) + 2ir/ A and 
A 2 = n\/ (cu>i) + n 3 j (clu 3 ) — n 4 / (0^4) + 27r/A are the 
phase mismatches for the lu 3 and LU4 SFG processes. 

Undepleted pump approximation. The coupled 
nonlinear equations (1) are often linearized assuming 
that the incident pump field E\ is much stronger than 
the other fields and therefore its amplitude is nearly con- 
stant (undepleted) during the evolution. Then Eqs. (1) 
are reduced to a system of three linear equations, 



id z A(z) = M(z)A(z), M = 



-Ai n; 







A 2 



(2) 



with Q p = Exy/U^, tt s = Exy/U^U, A(z) = 
[A2(z), A 3 (z),A 4 (z)] T , A 2 = g 1 ^ 2 e' Al V»3^ 4/2, A 3 = 
£1 £3-^2^4/2, A 4 = E 1 E 4 e- lA ^y / n 2 n 3 /2. Upon the 
substitution z — > t, Eq. (2) becomes identical to the time- 
dependent Schrodingcr equation for a three-state quan- 
tum system in the rotating-wave approximation, which 
is studied in great detail [8]; the vector A(z) and the 
driving matrix M correspond to the quantum state vec- 
tor and the Hamiltonian, respectively. We note that the 
quantity |A(z)| 2 = \A 2 (z)\ 2 + \A 3 (z)\ 2 + |A 4 (z)| 2 is con- 
served, like the total population in a coherently driven 
quantum system. By definition, in the adiabatic regime 
the system stays in an eigenvector of the "Hamiltonian" 
M. We assume that Ai(z) and A 2 (z) change linearly 
along z, which can be achieved, for example, by vary- 
ing A(z). Explicitly, we assume that either Ai = <5 — a 2 z 
and A2 = —5 — a 2 z, which is called "intuitive sweep" (for 
reasons that will becomes clear shortly) or Ai = 5 + a 2 z 
and A2 = — S + a 2 z which is called "counterintuitive 
sweep" . For the sake of generality, we take hereafter a 
as the unit of coupling and 1/a as the unit of length. 
Then the three eigenvalues of M will cross each other at 
three different distances z m (m = 1, 2, 3), thereby creat- 
ing a triangle crossing pattern [16-18]. These crossings 
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Fig. 1. (Color online) Sequential SFG of 0J4 field. Top 
frames: Diagonal elements (solid lines) and eigenval- 
ues (dashed lines) of the driving matrix M of Eq. (2) 
for the "intuitive" (left frames) and "counterintuitive" 
(right frames) phase mismatch sweep. The field intensi- 
ties are calculated numerically from Eqs. (1) for 6 = 2a, 
fii = O2 = ^3 = SI4 = a. Middle frames: undepleted 
pump, |£i(zi)| 2 = 100|£: 2 (z i )| 2 , with z { = -20a" 1 ; bot- 
tom frames: depleted pump, \Ei(zi)\ 2 = 2\E 2 (zi)\ 2 . 



allow us to design recipes for efficient broadband SFG, 
in analogy to adiabatic passage techniques in quantum 
physics [8,9,12,16-18]. Because of the analogy to the 
Schrodingcr equation the condition for adiabatic evolu- 
tion can be derived using the Landau-Zcner-Majorana 
model [19-21] and reads (for linear chirping and constant 
couplings): If^l > a, where Q x is the relevant coupling 
at the respective crossing. 

Figure 1 plots the eigenvalues of M of Eq. (2) vs z. 
Initially only the lo 2 field is present, hence the vector 
A = [A 2 , 0, 0]. If the evolution is adiabatic then there are 
two possible paths that the system can follow (marked 
by arrows) . If the phase match for the u 3 generation pro- 
cess occurs first (left frames of Fig. 1), then the energy 
is converted first to the u> 3 field and then to the W4 field. 
This "intuitive" two-step scheme extends the single-step 
adiabatic passage scenario for SFG [4,5]. Interestingly, 
we find that efficient energy transfer directly to the W4 
field is also possible through the "counterintuitive" di- 
rection of the local modulation period sweep when the 
phase match for the uj^ generation process occurs first 
(right frames of Fig. 1). Then the energy flows from the 
LU2 field to the lu^ field with almost no energy transferred 
to the intermediate co 3 field. 
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Fig. 2. (Color online) Efficiency of SFG of a; 4 field vs the 
couplings fii and fl 2 obtained by numerical integration 
of Eqs. (1) for "counterintuitive sweep" with S = a and 
O3 = O4 = a. Left frame: undepleted pump, |i?i(zi)| 2 = 
100|i?2(z 4 )| 2 , with z t = —20a" 1 ; right frame: depleted 
pump, lE^l 2 =2\E 2 (z t )\ 2 . 

Depleted pump. We have found by numerical inte- 
gration of the nonlinear system (1) that the described 
scheme is also applicable beyond the undepleted pump 
approximation, when the uj\ and lu 2 fields have compara- 
ble energies; this is demonstrated in the bottom frames of 
Fig. 1. Unfortunately, many optical parametric processes 
such as \bJ\ — lo 2 \, 2wi, 2lo 2 , wi + 2u;2 become possible 
in this case and it is not easy to find the conditions for 
broadband SFG of the oja field. 

The contour plot in Fig. 2 demonstrates the robustness 
of SFG of the o>4 field against parameter variations. SFG 
for an undepleted pump (left frame) is remarkably robust 
in confirmation of the simple analytic theory described 
above. SFG for a depleted pump (right) is less robust 
although relatively high SFG efficiency is still possible; 
because then the simple eigenvalue arguments cannot be 
used the interpretation is more difficult. 

Third harmonic generation. Third harmonic gen- 
eration is an important special case of SFG, which is 
readily treated in the adiabatic regime. The respective 
equations are derived from Eqs. (1), 
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= x (2) w/4cn w , fl 2uj = x {2) u\/3/4c^n 2uJ n 3u 
A u = E 1 ^/2n UJ /n 2ul , A 2ul = E 3 , A 3uj -- 



E 4 ^2n 3uj /3n 2uJ . 

Figure 3 shows numerical simulation of third harmonic 
generation. There are again two possible scenarios. If the 
phase match for the second harmonic generation process 
occurs first ("intuitive" sweep, left frame of Fig. 3), then 
the efficiency of the third harmonic is good, but we have 
some unwanted second harmonic left. For the "counter- 
intuitive" direction of the local modulation period sweep, 
when the phase match for the third harmonic generation 
occurs first, the second harmonic is strongly suppressed 
and a nearly complete transfer of energy to the third 
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Fig. 3. (Color online) Third harmonic generation for 
"intuitive" (left frames) and "counterintuitive" (right 
frames) phase mismatch sweep calculated numerically 
from Eqs. (1) for 6 = 4a, O w = 2a, £7 2tJ = 5a. 
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Fig. 4. (Color online) Efficiency of third harmonic gener- 
ation efficiency vs the couplings f2 w and f^w for "coun- 
terintuitive" sweep obtained by numerical integration of 
Eqs. (1) for S = 2a. 

harmonic takes place (right frame of Fig. 3). Figure 4 
demonstrates the robustness of the third harmonic gen- 
eration for the "counterintuitive" phase mismatch sweep, 
which indicates that this technique is broadband. 

Conclusion. We have used the analogy between the 
time-dependent Schrodinger equation and the SFG equa- 
tions in the undepleted pump approximation to propose 
an efficient broadband SFG technique. A local modu- 
lation period sweep along the light propagation creates 
crossings in the phase matching between different para- 
metric processes, which in combination with adiabatic 
evolution conditions allow efficient and robust SFG of 
the desired frequency W4 = 2wi +UJ2- While the physical 
picture is transparent in the undepleted pump approx- 
imation, the basic feature of the SFG process remain 
largely intact in the general regime of depleted pump. 
Chirped QPM gratings offer robustness against varia- 
tions of the parameters of both the crystal and the elec- 
tric fields, which include the crystal temperature, the 
wavelengths of the input electric fields, the crystal length 



and the angle of incidence. 

The present work can be viewed as a generalization of 
the idea of Suchowski et al. [4] from a single SFG to si- 
multaneous SFG processes in and beyond the undepleted 
pump approximation. This work is also a broadband al- 
ternative to the (narrowband) STIRAP-based third har- 
monic generation proposal of Longhi [3] . 

This work is supported by the European network 
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